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Rotating-wave approximation and its validity in multi-state quantum systems are studied through 
analytic approach. Their applicability is also verified from the viewpoint of generic states by the use 
of direct numerical integrations of the Schrodinger equation. First, we introduce an extension of the 
rotating-wave approximation for multi-state systems. Under an assumption that a smooth transition 
is induced by the optimal field, we obtain three types of analytic control fields and demonstrate 
their validity and deficiency for generic systems represented by random matrices. Through the 
comparison, we conclude that the analytic field based on our coarse-grained approach outperforms 
the other ones for generic quantum systems with a large number of states. Finally, the further 
extension of the analytic field is introduced for realistic chaotic systems and its validity is shown in 
banded random matrix systems. 



I. INTRODUCTION 

Theoretical and experimental studies of controlling 
quantum states have been attracting much attention be- 
cause of the theoretical progress in the field of quantum 
computing [l[ and of the technical developments in ma- 
nipulating atomic and molecular systems Q. Various 
control schemes are known in these fields: A 7r pulse 
in a transition between two eigenstates 3|, the nonadia- 
batic transitions induced by laser fields 4( , the STIRAP 
scheme by a counterintuitive pulse sequence for more 
than three- level systems etc. If we use an electron- 
ically excited state in the controlled system, a simple 
control of a pulse-timing can selectively break a chemical 
bond by pump and dump pulses @ . When two pathways 
exists from an initial state to a target state, quantum me- 
chanical interference can be utilized to modify the ratio 
of products, which is called a coherent control scheme 

These control schemes are effective for a certain class 
of processes but are not versatile for general multilevel- 
multilevel transitions. When we consider to design quan- 
tum devices with a large number of states interacting 
with a complex environment or with short-time laser 
pulses, such a multi-state control problem should be con- 
sidered. Optimal control theory [1, Q and genetic al- 
gorithms dl Gl| are most successful methods to solve 
this kind of complicated problems, but its implementa- 
tion and interpretation can be still difficult: What kinds 
of dynamical processes are involved in the controlled dy- 
namics? 

On the other hand, we know that highly excited sys- 
tems exhibit quantum chaotic features [121] . Such a 
" complex" quantum system driven by an external field 
[l3l Ibij is modeled by fully random matrix systems with 
a parameter, where statistical properties of eigenener- 
gies and eigenvectors [H, [la, [lj, LL8| are characterized 
under universality classes These statistical proper- 
ties stem from multilevel-multilevel interactions of eigen- 



states, which are related to the existence of many avoided 
crossings [2(1 [2l| ■ Hence it is necessary to consider the 
interaction between many eigenstates when we study dy- 
namics in such a system. G ong and Brumer applied the 
coherent control method 0, l22l ] to a quantum chaos sys- 
tem [23|, [24| and its prediction has been recently con- 
firmed by experiment. 

Several attempts have also been done to obtain con- 
trol fields analytically in multi-level systems [25], [26j] . We 
also derived an analytic optimal field to control the fully 
random matrix systems [2jj , which is based on the naive 
idea of "coarse-grained approach" [28|, |29[, and the re- 
sults are promising. In this paper, at first, we intro- 
duce several analytic schemes from the viewpoint of ap- 
plicability to generic quantum systems, and confirm the 
so-called rotating-wave approximation for multistate sys- 
tems. Next, we improve our previous approach [27| to 
deal with more realistic quantum systems with a banded 
random Hamiltonian, since the most realistic quantum 
systems may be modeled with such a banded random 
matrix [l9[ . The validity of the extended analytic field is 
evaluated by numerically solving the Schrodinger equa- 
tion for the multilevel-multilevel control problem. 



II. ANALYTIC FIELDS FOR MULTI-STATE 
SYSTEMS 



In a simple two-state transition problem, the Bloch 
vector in a three dimensional space is introduced to rep- 
resent the transition dynamics 3] . The rotating-wave ap- 
proximation (RWA) is also introduced if we note slowly 
varying dynamics on a rotating frame in the Bloch space, 
where a transition by a 7r-pulse is also represented as a 
rotation with an angle tt. The first problem considered 
in this section is "what is the natural extension of this 
representation in multi-state quantum systems." 
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A. The rotating-wave approximation 

Actual procedure of the RWA for two-state problems is 
to ignore off-resonant terms in the Schrodiner equation, 
where an intuitive interpretation of this is that a rapid 
motion induced by oscillating terms can be averaged out 
and only slow dynamics by near-resonant terms remain. 
We can extend this idea to the multi-state dynamics in- 
duced by an external field. 

We consider the Hamiltonian system 



H(t) = H + e(t)V 



(1) 



with an interaction term e(t)V by an external field s(t). 
If we introduce an eigenstate representation 



ih 



(2) 



the time evolution of {aj(t)} is represented by 

(Ej — E k )t 



a k (t) = ^ V kj<*j(t) ex P 



ift 



(3) 



where Ej and \ipj) are the j-th eigenvalue and eigenstate, 
respectively, and Vkj is a shorthand of (ip k \V\ipj). 

We consider the case [26| that the external field is a 
sum of oscillating terms with frequencies of the transition 
energies v jk = (Ej - E k )/h, 



3-k 



+ C.C, UJ jk 



E< - E k 



(4) 



According to the standard procedure of the rotating-wave 
approximation, we can separate the oscillating terms in 
the Schrodinger equation into slowly varying ones and 
rapidly changing ones. Further, we can introduce an 
approximation that all the oscillating terms are ignored 
when we consider the limit of infinitely long transition 
time under a non-degenerate condition between energy 
differences, 

Ej - E k ^ E y - E k , for j ^ j' and k ^ k' . (5) 



Note that the validity of this simplification of the 
Schrodinger equation is not trivial although this seems 
a direct extension of the usual RWA. 



B. Slow transition dynamics 

Consider a control problem from an initial state |$ ) 
at t = to a target state |$t) at t = T. The interaction 
representation [30( is introduced by 

|<M*)> - Ub(*. 0)|*o>, Ixo W) = U (t,T)\* T ), (6) 



where Uo(t2,ti) is a propagator by Hq from t = t% to 
t = t 2 . 

An overlap between |0o(i)} and |xo(i)} is parameter- 
ized by an angle (0 < < tt/2) with a phase a 
(0 < a < 2?r), 



(<M*)M*)> = ($o|(7o(0,T)|$t) =ie la sinO (7) 

which includes an orthogonal case (0 = 0). For the case 
of 0, the phase a is uniquely determined. Then, we 
can introduce orthogonal basis states, 



u aw u aw ic aw Ixoft)) - ie M sin0|</> o (t)) 

\<hW) = \Mt)), Xo(t)) = . 

COS0 

(8) 

From the result of optimally controlled dynamics in 
random matrix systems 27J, the time-evolution is as- 
sumed to be 



|V(t)> = \M*)) cos^ -te- ia \xo(t)} sin0, (9) 



in the limit of T — > oo. Substituting it into the 
Schrodinger equation, we obtain 



-\Mt))sm9~te- ia \x (t))cos9 



ih 



\4>o(t))cos9-ie- ia \ X o(t))smt 



(10) 



r 



If we operate (<j>o(t)\ from the left, the relation 



£(t) 

if) 



(Mt)\V\Mt)) cot 9 - ie- ia (j> (t)\V\ X o(t)) 



relation 
ih 



ie* a ( Xo (t)\V\4>o(t)) + (xoW|^|XoW)tan( 



(12) 



(11) It is almost impossible to satisfy these equations strictly 
is obtained, and the operation of (xo(^) gives another since such a field realizing the given dynamics Eq. §§§ does 
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not exist always. However, approximate fields can be 
obtained within a restriction to the form of Eq.(|3]). In 
the following, we try to determine e(t) under the RWA 
for several special cases. 



1 . Exact field for a transition from an eigenstate 



If we restrict e(t) to a sum of terms with transition fre- 
quencies from the 0-th eigenstate to the j-th eigenstate, 



e(t) = ^2s j0 e -^o(*- T ) + c.c, 



(21) 



The simplest example [26] is given by the case that 
$o) is an eigenstate, and |$t) is a linear combination of 
multiple eigenstates, 

|$o> = \<Pq), |$t) = sin @\<Po) +0088^^1^). (13) 
The angle 9 satisfies < 8 < 7r/2, and the normalization 

Ek-i 2 = 1 ( 14 ) 



we obtain the optimal field 



should be satisfied. We introduce 



(15) 



\Mt)) = W )e Eat/l \ 
|xo(t)) = sine|^o)e £o(t - T)/l?l 

+ coseErf J |^)e^( t - T ^, (16) 

and, from the relation Eq.®, the phase is determined by 

E T ir 
a= — -r 

Then, the orthogonal basis states are 



(17) 



\Mt)) = \w)e Eat/i \ |xo(t)> =J2 d ^)e EAt ~ T)/m - 

(18) 

The multiplication of (<y9o| from the left to Eq.fTU]) gives 
a relation 



.£(*) 
ih 



Vqq COt + Vojrfj 



-itv j0 {t-T) 



, (19) 



and the operation of (<£y| (j ^ 0) from the left gives 
another relation 



fc 

(20) 



(22) 



by ignoring all the oscillating terms (for details, see Ap- 
pendix^), where Q (> 0) is a constant or a slowly vary- 
ing function of time. If we substitute e e (t) into Eq. (fTO|) 
or Eq. ([20]) . 6 = f2 is easily shown. Since the target state 
|$t) is realized by 9 = | at t = T in Eq.©, the smallest 
rotation angle is § — 9. Thus, the smallest is deter- 
mined by 



n 



T. 



(23) 



which induces the perfect control (^tIV'W}! = 1 at the 
target time t = T. 

We note that e c (t) is obtained under the assumption 
of a slow transition Eq. ([9]) with a sufficiently long target 
time T. Then, the quantum state \ip{t)) driven by e e (t) 
stays in a plane determined by the two states, |0o(O) 
and Ixo(^)}? during the controlled dynamics from t = 
to t = T. 



2. Approximate field for a transition between multi-level 
states 



We study the case that the initial state |$o) and the 
target state are two different linear combinations of 
many eigenstates 

l*o) = 5>|^), l*r)=X)dilv J ->. (24) 



If we use Eq.® for the phase, the pair of orthogonal 
states is defined by 



\Mt)) = ^Cjl^e 



Ejt/ih 



Ixo(i)) - £ 



e - E * T/ih -ie ia Cj sinQ, 
cos 9 



Ejt/ih 



(25) 
(26) 
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Based on the assumption that the controlled state rep- 
resents a smooth rotation Eq.© and that e(f) has a re- 
stricted form Eq.Q, we can define 



ff a (t) = nn 

j,k (#i) 



C J a k c ia c -iuj ik t 



V, 



jk 



(27) 



as an approximate control field (for details, see Appendix 
HU). If we substitute e a {t) into Eq.JIl]) and Eq.(T2]), we 
obtain 



9 = fl[l + 0(N- 



(28) 



by ignoring all the oscillating terms, where N is the num- 
ber of eigenstates contained in the state, and {c, } and 
{dk} are assumed as random complex numbers. Note 
that the approximate field s a (t) is valid in the limit of 
N — » oo as well as T — > oo. 



3. Another field for a transition containing many 
eigenstates 

We use another approach to the approximate analytic 
field in the limit of infinitely many eigenstates. By the 
assistance of the optimal control theory [271 ] , we can in- 
troduce a control field defined by 



where 



. , 2?ifi 
ScsW — ,Re 



'(^o(*)|V|xo(t)> 



(0o(*)l^ 2 |£o(t)) 



(It. 



(29) 



(30) 



Substituting e cg (t) into Eqs. (fTTjl and (fT2|) . we obtain 



6{T) = / 0(i)df w riT 
Jo 



(31) 



under the limit of T — > oo and 2V — > oo. This is the 
optimal field based on the coarse-grained approach we 
derived before [Si Hi]. 



III. NUMERICAL EVALUATION FOR 
GENERIC SYSTEMS 

In the previous section, we have obtained several op- 
timal fields for multi-state transitions, i.e., e c (t), e a (t), 
and £ cg (t). Applicability of these fields depends mainly 
on the validity of the RWA. Furthermore, for the cases 
of e e (t) and £ a (£), the matrix element Vjk in the denom- 
inator may deteriorate the results while the effect of the 
number of states should be checked for e a (t) and s cg (t). 
All these points can be verified in a direct numerical inte- 
gration of the Schrodinger equation for generic quantum 



systems with a random matrix Hamiltonian under Gaus- 
sian Orthogonal Ensemble (GOE) or Gaussian Unitary 
Ensemble (GUE) [H. 

The GOE random matrix represents the case of "com- 
plex'' Hamiltonian systems with a time-reversal symme- 
try [li]. This is constructed by a real symmetric ma- 
trix with off-diagonal elements v subject to a distribution 
function 



(v)dv oc exp 



2A 2 



dv. 



(32) 



This means that the density of v has a maximum value 
at v = 0. 

The GUE random matrix represents Hamiltonian sys- 
tems without time-reversal symmetry [19J. This is con- 
structed by a Hermitian matrix with off-diagonal ele- 
ments z = x + iy subject to a distribution function 



Pgue(z)dz oc exp 



2A 2 



dz 



(33) 



If we take a polar-angle representation (r, (f>), z = x+iy = 
re 1 ^ , the distribution function for r is 



P sue (r)dr oc r exp 



2A 2 



dr 



(34) 



This means that the density of r vanishes at r = 0. Thus, 
the number of small off-diagonal elements in GUE ma- 
trices are relatively fewer than that in GOE matrices. 



A. Final overlaps of transition from an eigenstate 

For the numerical evaluation, Hamiltonian matrices 
Hq and V are created by random numbers subject to 
GOE or GUE. Hq is scaled so that the average spacing 
AE = (\Ej — Ej-\\) is unity, and the target time T is 
shown in units of h/ AE. The exact control field e c (t) 
is defined immediately after the target state |$t), i-e., 
coefficients {dj}, is determined. Then, the final overlap 
|($t|V'(T 1 ))I i s obtained by numerical integration of the 
controlled dynamics by e e (t) written in the Schrodinger 
equation without RWA. FigfT] represents the result ob- 
tained as an ensemble average over 100 different samples 
of Ho, V, and {dj}. 

It is shown that e c (t) works only for sufficiently large T 
in spite of the exact control field under the RWA, which 
means that the RWA is valid only for such a large T. 
The reason why the RWA breaks down for smaller T is 
that the field amplitude for the frequency corresponding 
to the — > j transition becomes accidentally large when 
V}o| ~ 0. In order to avoid the breakdown, we must 
use the field e e (t) with the smaller amplitude which is 
realized by a larger T. If we take the limit T — » oo, in 
this case, it is expected that \{^t\^{T)) \ — > 1 since e (t) 
is the exact control field within the RWA. 
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FIG. 1: Final overlaps by e c (t) for the transition from an 
eigenstate to generic states constructed from 5 (dashed curve), 
10 (thin curve), and 20 eigenstates (thick curve) in (a) GOE 
and (b) GUE random matrix systems. Each curve is obtained 
as an average over 100 different random configurations of Hq, 
V, and {dj}. 



FIG. 2: Final overlaps by £&(t) for the transition between 
generic states constructed from 5 (dashed curve), 10 (thin 
curve), and 20 eigenstates (thick curve) in (a) GOE and (b) 
GUE random matrix systems. Each curve is obtained as an 
average over 100 different random configurations of Ho, V, 
{cj}, and {dj}. 



Difference in the performance of this field with respect 
to the universality class is also explained by the differ- 
ence in probability distribution functions for smaller \Vjo\ 
(See Eq. ([32|) and Eq. ([34f ). When we use larger matrices, 
larger target times will be necessary to obtain a certain fi- 
nal overlap since the probability for accidental divergence 
of the amplitude tends to occur for many off-diagonal el- 
ements. 



B. Final overlaps in transition between multi- level 
states 

In case of multilevel-multilevel transitions, the same 
procedure can be executed with a numerical integra- 
tion of the Schrodinger equation without the RWA. In 
this case, control fields e a (t) and s cg (t) is determined by 
Hq, V, {cj}, and {dj}, which are created from pseudo- 
random numbers according to distribution functions un- 
der each universality class. Fig [2] shows the performance 
of e a (t) for GOE and GUE random matrices. Almost 
the same properties can be observed as the case of the 
transition from an eigenstate. However, even if we take 
the limit T — ► oo, the final overlap |( ( I ) t|V'(^ 1 )}| does not 
seem to converge to 1. This is because e a (t) is an ap- 
proximate field obtained by ignoring quantities with the 



order of 1/N. 

Figgis the result by £ cg (t). In GOE case (FigG^a)), 
much improvement can be seen when it is compared 
to one by e a (t), while it is almost the same for GUE 
(FigEHb)). The overlap does not converge to 1 in the 
limit of T — * oo since s cg (t) is valid under N — > oo. The 
important thing observed here is that, in case of e cg (t), 
longer target times are not necessary even for larger ma- 
trices. This is significant in the practical application of 
controlling large systems. 



IV. EXTENSION TO REALISTIC CHAOTIC 
SYSTEMS 

The full random matrix systems studied in the previous 
section is a model of strongly chaotic systems. In order 
to consider realistic quantum systems, one of candidates 
to be studied is a banded random matrix [3l| . In this 
section, we try to improve our external field e cg (t) for 
the case of the banded random matrix 13211 . 
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FIG. 3: Final overlaps by e cs (t) are shown for the transition 
between generic states constructed from 20 (thick curve), 50 
(thin curve), and 100 eigenstates (dashed curve) in (a) GOE 
and (b) GUE random matrix systems. Each curve is obtained 
as an average over 100 different random configurations of Hq, 
V, {cj}, and {dj}. 
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FIG. 4: Examples of the interaction Hamiltonian V (50 x 
50). The magnitude of matrix elements are shown for (a) full 
random matrix subject to GOE, and (b) a banded random 
matrix subject to GUE (Ao = 15). 



A. Analytic field for banded random matrix 
systems 



We consider the case that the interaction Hamiltonian 
V is a banded random matrix in the eigenstate represen- 
tation of Hq. The elements of V are random complex 



numbers with the following distribution 



V S k\' 



(h\v\<p k 



exp 



A 2 



We introduce a new field £brm(i) as an extension of £ C g(^)> 

£br,„(t) = Re [ A jkC*V jk d k e lw ^] , (36) 
jk 

with an extra-amplitude factor Aj k . The coefficients 
cij(t) in Eq.© satisfy the Schrodinger equation 

ifi-^a k (t) = [A]kc*d k + A* k] c k d*} \V jk \ 2 aj(t) 

after ignoring all the oscillating terms. If we assume that 
the transition is smooth, aj (t) should be written as 



a k (t) = c k cos 



7rf \ 

2T J 



id k sin 



TTt \ 

2TJ ' 



(38) 



where we take O = in Eq.Q for simplicity. Substitut- 
ing these coefficients into Eq. ([3T|) . we obtain a relation 



iirh 
2T 



"' 7, ' shl ' if) ~ idkcos (if 



It, 



Ajk\cj\ d k cos 



TTt 

2T 



iA* kj c k Irfj T sin ( ] 



2T I 



Jk 



(39) 



under the assumption of random phases [27] . Finally, we 
obtain the conditions for the coefficients Aj k 



Ajk = A M ~ — exp 



T 



(Ej — E k ) 2 



A Z 



(40) 



If we consider the case that those coefficients Cj and d k 
have Gaussian distribution functions in the energy space, 



(N 2 ) oc exp 
(\d k \ 2 ) oc exp 



(E j - E c f 



(E k - E d f 



(41) 



with centers E c and Ed and widths A c and Ad, we can 
define 



itTi 



£brm(t) = y^Re [c*V jk d k e"^ 4 ] exp 



jk 



(Ej - E k ) 2 



This field has a finite amplitude only when 
A c < A and A d < A . 



(42) 



(43) 



If not, the field has an infinite amplitude in the limit 
of Ej,E k — > ±00 due to the exponential factor Aj k in 
Eq. ([42|) . Thus, the analytic field is refined when the 
widths of the initial and target states are relatively small 
compared to the width of the banded random matrix el- 
ements. 
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FIG. 5: Improvements by Ebrm(i) is shown for (a) GOE case 
and (b) GUE case. The initial and the target states are both 
Gaussian wavepacket in the energy space with A c = A d = 15. 



B. Numerical evaluation 



FIG. 6: Improved final overlap by Ebrm(i) (thick curve) is 
compared to that by e cg (t) (dashed curve) for the case of 
banded random matrix systems with A c = Ad = Ao = 15 
subject to (a) GOE and (b) GUE. The result by e cg (t) in a 
full random matrix system is also shown for comparison (thin 
curve). Each curve is obtained as an average over 100 different 
random configurations of Ho, V, {cj}, and {dj}. 



We confirm the validity of £brm (t) for the system with a 
banded random-matrix interaction. The numerical test 
is configured as follows. The initial and target states 
are defined as quantum vectors with random complex 
coefficients Cj and dj subject to Eq. ljlTj) . Here, we choose 
A c = A d = 15, E c = -10, and E d = 10, where H 
is a 50 x 50 random matrix of GOE and GUE, and is 
scaled so that its eigenvalues {Ej} are distributed in an 
interval [—25, 25]. The interaction Hamiltonian V is also 
a 50 x 50 matrix while its elements obey a banded-random 
distribution in the eigenstate representation of Ho with 
A . 

The optimal field Eq. (|42|) is calculated from those 
quantities {cj}, {dk}, {Vjk}, and {Ej} with parame- 
ters T and Ao- In order to check the validity of £t,rm(i) 
Eq. (|42p . we solve the initial value problems with Hamil- 
tonian Eq.([T]) driven by Eq. (|42[) for various band widths 
Ao of the interaction Hamiltonian V. The results are 
shown in Figure [5] When we use the original analytic 
field s cg (t) Eq. (|2T)]) . the performance of the optimal field 
(dashed curve) decreases for the banded matrices with 
smaller widths. On the other hand, the final overlaps 
(solid curve) by the refined analytic field Eq. (j42|) does 
not change even for the smaller width untill the limit 
Ao w A c or Ao ~ A^. In Fig[5J the performance of 
£brm(t) is shown for the various target time T . The final 



overlap by £brm(i) (thick curve) is much improved from 
the result by s C g(t) (dashed curve), which is comparable 
to the level of the original result by a full random matrix 
(thin curve). 

V. CONCLUSION 

We studied several forms of analytic fields to steer 
quantum states in generic systems. Through numerical 
evaluations of those fields, the rotating-wave approxima- 
tion for multi-level systems are also validated in the limit 
of the large target time. It was shown that the ana- 
lytic field defined through the coarse-grained idea can be 
used in the limit of many states. We extended our pre- 
vious analytic approaches for controlling complex quan- 
tum systems to deal with more realistic systems with a 
banded random matrix. The key ingredient is the ampli- 
tude factor Ajk , which is an exponentially growing func- 
tion, introduced in the analytic optimal field Eq. (|4"2")l . We 
showed that the new analytic field outperforms the pre- 
vious field Eq. (|29p for a full random matrix. In the near 
future we will apply this optimal field to quantum chaos 
systems such as quantum kicked rotors (top) [33j and to 
more realistic molecular systems [34], [35| . 
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APPENDIX A: DERIVATION OF THE EXACT FIELD e e (t) 



Substitution of the restricted form Eq. ([2"Tj) into Eqs. (fT9|) and (|20|) gives relations 



in ^ V 



Vqq cot 9 + Y Vojdj e ~^°( t - T ) 
j¥0 



(Al) 



and 



6» <L 



(A2) 



respectively. By dropping all the oscillating terms from 
these expressions, we obtain 

ih 9 = — Vojdj £* , ift. d/fl = Vjo Ejo (3) 

as relations under RWA. From the second relation, Sjo is 
determined to be 

e , = g-rxi, (4) 

and it is easily shown that this definition also satisfies 
the first relation if we note Voj — V* Q and \dj\ 2 = 1- 
Thus, the exact control field under RWA is obtained as 
Eq.£B). 



II. DERIVATION OF THE APPROXIMATE 
FIELD £a(t) 

By substituting Eqs.([2l)|) and (US]) into Eas.fTT]) and 
(fT^)l . we obtain equations for e(t) in the eigenstate rep- 
resentation, 

9 = - £ -^J2 c j{ CkCot9 - ie ~' a ^) V Jk eiUJlkt > (!) 

= ^ E J J ( ie< ° c * + tan 0) y jfe e*""*, (2) 

respectively. If we assume the restricted form Eq.(|4]) as 
the control field, these equations can be simplified to 

° = --k {e ]k +el j )c*(c k cot9-ie~ ia d k )v ]k , 

(3) 

Q=^Y. {sjk + e* kj )d*^e ia c k +d k t a ne)v jk , 

(4) 



after dropping all the oscillating terms. The problem is to 
determine the amplitude factors Ej k which satisfy these 
equations for given 9. 

Here, we try to define a simple solution that is inde- 
pendent of 9. Such a solution, if exist, should satisfies 
Eqs.© and (U|) for any value of 9. From Eq.© with 
9 = 7r/2 or Eq.((4|) with 9 = 0, a condition for Sj k , 

d =-f- X, {ejk + £* kj ) c*jV jk d k (5) 

is required. While it is still difficult to define £j k satisfy- 
ing this equation exactly, an approximate one is defined 
when the number of states N is large enough. By intro- 
ducing amplitudes, 

ia 3* 

e 3k = e -^nn, (6) 

Vjk 

as an extension of Eq. (01, the required Eq.© leads to 
Eq.(28J by the use of 

£ c* Cj d* k d k = 1 + OiN- 1 ), 

'Vlj'vh.. = OiN- 1 ), (7) 

iMft) 

under the assumption that there are no special correla- 
tions between Cj and dj other than the orthogonal rela- 
tion J2j Cjdj = 0. 

Moreover, Eq.© approximately satisfies Eqs.Q and 
((4]) for any 9, because small quantities 

Y c*jCjd* k c k , Y d*d jC * k d k , (8) 

of the order 0(N~ 1 ) are also ignored in the limit of N — > 
00. Thus, the control field Eg. ([27)1 is approximately valid 
for a transition between multi-level states. 
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